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We report measurements of CP violation parameters in — > K^n'^'y transitions based on a data 
sample of 535 x lO^B-B pairs collected with the Belle detector at the KEKB asymmetric-energy 
e'^e~ collider. One neutral B meson is fully reconstructed in the K%Ti^'y mode. The flavor 

of the accompanying B meson is identified from its decay products. We obtain time-dependent 
and direct CP violation parameters S and A for a JfgTr" invariant mass up to 1.8GeV/c^ as 
Sk%^o-, = -0.10 ± 0.31 ± 0.07 and ^^g^o^ = "0.20 ± 0.20 ± 0.06. For a K%-k° invariant mass near 
the A"*" (892) resonance, we obtain Sk'O-, = -0.32lo ± 0.05 and Ak'O^ = -0.20 ± 0.24 ± 0.05. 

PACS numbers: 11.30.Er, 13.25.Hw 



The radiative b ^ sj penguin process is sensitive 
to physics beyond the standard model (SM), and time- 
dependent CP violation in decays of the type — > 
fcpl, where fcp is a CP eigenstate, has drawn much 
theoretical and experimental interest recently [l|, 0, H, 0j 
[i,[i,0,[3- Within the SM, the photon emitted from a B° 
{B^) meson is predominantly right-handed (left-handed). 
A flip of photon polarization is suppressed by the quark 
mass ratio 2ms/mb [H. Hence, for fcp = K*'^ (— > i^^Tr^) 
the SM predicts a small time-dependent CP asymmetry, 
which arises from the interference between decay ampli- 
tudes with and without B'^-B'^ mixing. The same sup- 
pression is expected for final states of the type i?" 
P^Q'^^ where P" and are any C eigenstate spin-0 
neutral particle (e.g. P" = and Q° = 7r°). However, 
there are estimates predicting an enhancement of the 
asymmetry up to 0.1 due to strong interactions For 
the case of P*' — > K*^ (— > K^n^)^, explicit computations 
support a small asymmetry, S = —(3.5 ± 1.7) x 10^^ 0| 
or5 = -(2.2±1.2j;? o)xlO-2 in the SM. A significant 
deviation from the small SM expectation could indicate 
new physics. 

Since the time-dependent CP asymmetry is not ex- 
pected to change significantly as a function of K^tt^ in- 
variant mass {Mj^o^^o ) (3] , we perform two measurements: 
one for P" K*^{-^ KI'k^)-^ Q by requiring M^^-o^o to 
lie in the range 0.8 GeV/c^ < M^o^o < l.OGeV/c^, and 
the other for the full range of M^o^o below 1.8GeV/c^. 
For simplicity, we refer to these two analyses as K*^^ and 
irg7r°7, respectively. The measurement for if *°7 is the- 
oretically cleaner than the measurement for if^ 71*^7, but 



the latter has more statistical power. Similar measure- 
ments have been previously reported by both Belle Q 
and BaBar (i| based on 275 x 10^ and 232 x 10^ PP 
pairs, respectively. In this Communication, we update 
the measurements of CP parameters for P° K*^^ 
and P° iCg7r°7 based on a data sample that contains 
535 X 10^ PP pairs. 

At the KEKB asymmetric-energy e^e~ (3.5 on 
8.0 GeV) collider 0, the T(45') is produced with a 
Lorentz boost of /?7 = 0.425 along the z axis, which 
is defined as the direction antiparallcl to the beam 
direction. In the decay chain T(45') B^B^ /icc/tag, 
where one of the P mesons decays at time tree to a final 
state /roc, which is our signal mode, and the other decays 
at time ttag to a final state /tag that distinguishes between 
P" and P°, the decay rate has a time dependence given 
by 
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S sin(Amd At) 



^cos(AmrfAt) 



(1) 



Here S and A are CP-violation parameters, tbo is the 
P° lifetime, lS.m^ is the mass difference between the two 
P° mass eigenstates, Ai is the time difference tree — ^tagj 
and the 6- flavor charge q = +1 (—1) when the tagging 
P meson is a P° (P°). Since the P° and P° mesons are 
approximately at rest in the T(4S') center-of-mass system 
(c.m.s.). At can be determined from the displacement in 
z between the /roc and /tag decay vertices: At ~ (zroc — 
ztag)/(/37c) = Az/(/37c). 
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The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber, an array of aerogel 
threshold Cerenkov counters, a barrel- like arrangement 
of time-of-flight scintillation counters, and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals lo- 
cated inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. An iron flux-return located out- 
side of the coil is instrumented to detect if ° mesons and 
to identify muons. The detector is described in detail 
elsewhere [llj. Two inner detector configurations were 
used. A 2.0 cm radius beampipe and a three layer silicon 
vertex detector (SVDl) were used for the first sample of 
152 X 10^ BB pairs, while a 1.5 cm radius beampipe, a 
four layer silicon detector (SVD2) and a small-cell inner 
drift chamber were used to record the remaining 383 x 10^ 
BB pairs [H. 

For high energy prompt photons, we select the cluster 
in the ECL with the highest energy in the c.m.s. from 
clusters that have no associated charged track. We re- 
quire 1.4 GeV < E^-'"^-^- < 3.4 GeV. For the selected 
photon, we also require Eg/E25 > 0.95, where Eg/E25 
is the ratio of energies summed in 3 x 3 and 5x5 ar- 
rays of CsI(Tl) crystals around the center of the shower. 
In order to reduce the background from tt*^ or ry decays, 
photons from candidate tt*^ ^ 77 or 77 ^ 77 decays 
are rejected using a likelihood described in detail else- 
where [l3|. The polar angle of the photon direction in 
the laboratory frame is restricted to the barrel region 



require that the 7r+ candidate is not positively identified 
as any other particle species {K~^, p+, e+ and /x"'"). In 
the i?° A'+7r~7 sample, the A'"*" candidate is selected 
from charged tracks identified as kaons, and the 7r~ can- 
didate from the rest of the tracks. 

The B" A:°7r°7 and B+ K"gTT+j modes are 
reconstructed simultaneously and a single candidate is 
selected from possible multiple candidates amongst the 
two modes in order to reduce the cross-feed background 



from 5+ — > ifg7r+7 in B^ ~> K'^tt'^j, where a tt^ is se- 
lected instead of the 7r+. The best candidate selection 
is based on a likelihood ratio TZ calculated from likeli- 
hood variables for signal {Csig) and background {Chkg) 
as 72. = Csigl {Csig +Cbkg), where the likelihood variables 
are obtained from a Fisher discriminant J-^ [ist . which 
uses the modified Fox- Wolfram moments [lij as discrim- 
inating variables. Hereafter, we denote the likelihood 
ratio with the likelihood variable name in parentheses, 
as TZ{^) in this case. We select the candidate that has 
the largest TZ{T). We form two kinematic variables: the 
energy difference AE 
energy constrained mass Mb 



- - Eg;!"^- and the beam- 

(pc^m.s.)2^ 

where E^-^^^- is the beam energy, and i^^™-^- and pg™-^- 
are the energy and the momentum of the candidate in the 
c.m.s. The signal region in AE and Mhc, which is used 
for the measurements of CP-violating parameters, is de- 
fined as -0.2 GeV < AE < 0.1 GeV and 5.27GeV/c2 < 
Afbc < 5.29GeV/c^. In order to determine the AE- 
A/bc dependent signal fraction, a larger fitting region, 

is 



of the ECL (33° < 9^ < 128°) for SVDl data, but is -0.3 GeV < AE < 0.5 GeV and 5.2GeV/c2 < Mtc 
extended to the end-cap regions (17° < 9j < 150°) for 
SVD2 data due to the reduced material in front of the 
ECL. 

Neutral kaons {Kg) are reconstructed from two oppo- 
sitely charged pions that have an invariant mass within 



±6MeV/c^ (2(t) of the mass [ij. The tt+tt" vertex 



is required to be displaced from the interaction point (IP) 
by a minimum transverse distance of 0.22 cm for high mo- 
mentum (> 1.5GeV/c) candidates and 0.08 cm for those 
with momentum less than 1.5GcV/c. The direction of 
the pion pair momentum must agree with the direction 
defined by the IP and the two pion vertex point within 
0.03 rad for high-momentum candidates, and within 0.1 
rad for the remaining candidates. Neutral pions {tt'^) arc 
formed from two photons with an invariant mass within 
±16MeV/c^ (3(t) of the tt" mass. The photon momenta 
are then recalculated with a tt" mass constraint. We then 
rec^uirc the momentum of tt" candidates in the c.m.s. to 
be greater than 0.3GeV/c. The K'^tt'^ invariant mass, 
M^o^o, is required to be less than 1.8GeV/c^. 

We also reconstruct _B" /v +7r^7 and i?+ /v^7r+7 
candidates as control samples in a similar way. Charged 
tracks other than those from Kg are required to origi- 
nate from the IP (within 5 cm in z and 1.4 cm in r-0); 
the transverse momentum {pt) is required to be greater 
than O.lGeV/c. In the i?+ KgTr~^j sample, we also 



used. 

In order to suppress the background contribution 
from continuum light quark pair production processes 
(e+e~ qq with q ~ u,d,s,c), which we here- 
after refer to as qq, we form another likelihood ratio 
TZ{T, cos 9b, cos Oh) by combining T with cos9b and 
cos 9h, where 9b is the polar angle of the B meson can- 
didate momentum in the c.m.s. and 9h is the helicity 
angle defined as the kaon momentum direction with re- 
spect to the opposite of the B momentum in the K-n 
rest frame. The helicity distributions for signal and back- 
ground are determined from the K^tt~^ sample 
for three mass regions: 0.8GeV/c^ < Mktt < l.OGeV/c^ 
(MRl), 1.3GeV/c2 < Mk^ < 1.55GeV/c2 (MR2), and 
the remaining range up to 1.8GeV/c^ (MR3). The dif- 
ference of background distributions in B^ — > K'^tt^'j and 
B^ KgTT^j decay modes is corrected using sideband 
data samples. In addition, a helicity dependent efficiency 
correction is apphed. Specific TZ{T, cos 9b, cos 9 h) selec- 
tion criteria are applied depending on both the mass re- 
gion and the flavor tagging information. Background 
contributions from B decays, which are considerably 
smaller than qq, are dominated by cross-feed from other 
radiative B decays. 

The 6-flavor of the accompanying B meson is identified 
from inclusive properties of particles that are not associ- 
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ated with the reconstructed signal decay. The algorithm 
for flavor tagging is described in detail elsewhere 17| . We 
use two parameters, q defined in Eq. ^ and r, to rep- 
resent the tagging information. The parameter r is an 
event-by-event flavor-tagging quality factor that ranges 
from to 1: r = when there is no flavor discrimina- 
tion and r = 1 when the flavor assignment is unambigu- 
ous. The value of r is determined by using Monte Carlo 
(MC) and is only used to sort data into seven r intervals. 
Events with r > 0.1 are sorted into six r intervals. The 
wrong tag fraction w and the difference Aw in w between 
the B° and decays are determined for each of the six 
r intervals from high-statistics control samples of semi- 
leptonic and hadronic b ^ c decays. If r is less than or 
equal to 0.1, we set w to 0.5, and therefore the accompa- 
nying B meson does not provide tagging information in 
this case. 

The vertex position of the signal-side decay of B'^ 
Kgir'^'y and 5+ A'^7r+7 is reconstructed from the Kg 
trajectory with a constraint on the IP; the IP profile 
(ax — 100 fim, (Tj^ ~ 5 fim) is smeared by the finite B 
flight length in the plane perpendicular to the z axis. 
The Kg vertex is displaced from the B vertex and often 
lies outside of the SVD. In this case, the vertex resolution 
is not good enough for a time-dependent CP asymmetry 
measurement. Therefore, both pions from the Kg decay 
are required to have enough hits in the SVD: at least 
one layer with hits on both the z and r-cj) sides and at 
least one additional hit in the z side of the other layers 
for SVDl, and at least two layers with hits on both sides 
for SVD2. The other (tag-side) B vertex is determined 
from well reconstructed tracks that are not assigned to 
the signal side. A constraint to the IP profile is also 
imposed. 

After all the selections are applied, we obtain 4078 can- 
didates in the Ai^-Mbc fit region, of which 406 are in the 
signal box. We perform an unbinncd maximum likelihood 
(UML) fit to the AE-Mbc distribution in order to re- 
solve signal, BB background and qq background compo- 
nents. The signal probability density function (PDF) is 
obtained from MC. We use a two-dimensional histogram 
of MC simulated data, for which the peak position and 
the width are corrected to account for differences between 
data and simulation using the 5° — > if+7r~7 control 
sample. The two-dimensional PDF for the BB back- 
ground, which populates more the lower AE region, is 
obtained from MC. For qq background, we use the prod- 
uct of two one-dimensional PDFs: the ARGUS param- 
eterization [l3| for Mbc and a second order polynomial 
for AE. Five parameters, which are the signal fraction, 
the BB background fraction, ARGUS shape parameter 
a, and two polynomial coefficients (ci, C2), are the free 
parameters in the fit. 

We first fit the entire M^^o^o region, and then fit the 
MRl, MR2 and MRS mass regions separately with the 
three background shape parameters (a, ci and C2) fixed 
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FIG. 1: (a) AE distribution within the Mbc signal slice and 
(b) Mbc distribution within the AE signal slice for the whole 
Mf^o^o region. Points with error bars are measured data. The 
solid curves show the fit results. The dotted curves show the 
qq background contributions, while the dashed curves show 
the sum of qq and BB background contributions. 



at the values obtained from the fit to the full range. 
Figure [T] shows a AE (Mbc) projection of the fit result 
in the Aibc (AE) signal slice for the entire Mj^o^o re- 
gion. For the MRl, MR2, and MR3 samples, we find 
112.5 ±12.0, 28.7 ±7.1, and 35.2 ± 10.0 signal events, 
respectively, with signal-to-background ratios (S/N) of 
1.91, 0.81, and 0.35. The level of BB background is 
only 7.7 ±4.2, 4.1 ±2.4, and 5.3 ± 5.2 BB events, re- 
spectively. As expected, in the K*^ mass region the S/N 
ratio is high. 

We determine S and A from an UML fit to the ob- 
served At distribution. The PDF expected for the sig- 
nal distribution, ■psig(Ai; 5, A, q, w, Aw), is given by the 
time-dependent decay rate [Eq. ([1])], modified to incor- 
porate the effect of incorrect flavor assignment; the pa- 
rameters Tro and Arrid are fixed at their world-average 
values 



14| . The distribution is then convolved with the 



proper-time interval resolution function i?sig , which takes 
into account the finite vertex resolution. The parameter- 
ization of i?sig is the same as that in the previous mea- 
surement [3l , while the parameter values are updated for 
the whole dataset. The PDF for BB background events 
{Vg-g) is modeled in the same way as signal, but with 
different lifetime and CP-violating parameters, while the 
resolution function Rgg is the same as i?sig. The effec- 
tive lifetime of BB background is obtained from a fit to 
the MC sample; the result is 1.356 ± 0.045 ps. The back- 
ground is assumed to have no CP asymmetry. Since 20% 
of the BB background are from non-radiative B^ decay, 
assigning the maximum asymmetries to this component, 
possible CP asymmetries in the background [S = ±0.2 
and A = ±0.2) are taken into account in the systematic 
error. 

The PDF for qq background events, Vqq, is modeled 
as a sum of exponential and prompt components, and 
is convolved with a double Gaussian which represents 
the resolution function Rqq. All parameters in Vqq and 
Rqq are determined by a fit to the At distribution of a 
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background-enhanced sample in the AE-M\^c sideband 
region. 

For each event, the following likelihood function is eval- 
uated: 



/ + 00 
-OO 



UgVsig{At')R,ig{AU ' At') 
fgj^Pgj^{At')Rgj^iAU - A^') 
(1 - /sig - fsB)P,^{At')Rgg{At, - At' 
/oiFoi(At,), 



d{At') 
(2) 



where Pqi is a Gaussian function that represents a small 
outlier component with fraction /oi [l^ . The probabil- 
ity of signal (/sig) and background {fgg) are calculated 
on an event-by-event basis using the results of the two- 
dimensional AE-M\^c fit, and are then multiplied by a 
factor that depends on the flavor tagging bin r. The r 
distributions of signal and qq background events are es- 
timated by repeating the AE-M^c fit procedure for each 
r interval with the three background shape parameters 
fixed to the full range result. The BB background dis- 
tribution is estimated from MC since the number of BB 
background events in data is limited. 

The only free parameters in the CP fit to B^ K^tt'-'j 
are S^o-^o^ and Ak'^.-^o^, which are determined by max- 
imizing the likelihood function L = J^^ Pj(Ati; 5, ^), 
where the product is over all events. We obtain 

Sro^o^ = -0.10±0.31(stat) ±0.07(syst), (3) 
■^K°^o^ = -0.20±0.20(stat) ±0.06(syst), (4) 

where the systematic errors are obtained as discussed be- 
low. We define the raw asymmetry in each At bin by 
{Ng=+i - Ng=^i)/{Ng=+i + iV,=_i), whcrc A^,^+i 
is the number of observed candidates with q = +1 (—1). 
Figure [2] shows the At distributions of the events with 
0.5 < r < 1.0 for q = +1 and q = —1 and the raw 
asymmetry. 

We perform the following fits to confirm the validity of 
our procedure: a i?+ lifetime fit for the B^ A'§7r+7 
sample gives 1.49 ± 0.10 ps, which is consistent with the 
nominal B^ lifetime [lj|; a B^ lifetime fit for the _B" 
Kgir'^j sample gives 1.53j^n i? ps, which is consistent with 
the nominal B^ lifetime [1J| ; and a CP asymmetry fit for 
the B^ A'g7r+7 sample gives an asymmetry consis- 
tent with zero {S = 0.20 ± 0.18, A = 0.05 ± 0.11). This 
is expected since the charged decay is completely flavor 
specific regardless of the photon polarization. A fit to 
MRl data gives ^^^.o..^ = -0.32;°J^(stat) ± 0.05(syst), 
Ak'O^ -0.20±0.24(stat) ±0.05(syst). A fit to non- 
MRl data gives S ^ +0.50 ± 0.61(stat) ± 0.29(syst) and 
A = -0.20±0.37(stat) ±0.13(syst). These results are 
consistent with those from the full Mj^o^o sample. 
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FIG. 2: (Top) Proper time distributions for B° K°sit°^ 
for g = +1 (left) and g = -1 (right) with 0.5 < r < 1.0. 
The solid curve shows the total and dashed curve shows the 
signal component. (Bottom) Asymmetry in each At bin with 
0.5 < r < 1.0. The solid curve shows the result of the UML 
fit. 



We evaluate systematic uncertainties in the following 
categories by fitting the data with each fixed parameter 
shifted by its error: uncertainties from physics parame- 
ters such as Arud, t^o, effective lifetime and CP asym- 
metry of BB background, uncertainties in the knowledge 
of qq background At PDF, uncertainties in the flavor 
tagging, uncertainties in the signal and background frac- 
tions, uncertainties in the resolution functions. A possi- 
ble bias in the flt is checked by fitting a large MC sample. 
The fit result is consistent with the input value within 
the statistical error. We quote this statistical error as 
the possible fit bias. Uncertainty due to the vertex re- 
construction is estimated by using the B^ J/ipKg 
control sample, where the tracks from J/ip are ignored 
for the study of the signal side vertex reconstruction us- 
ing the Kg trajectory. The effect of SVD misalignment 
is estimated by artificially displacing the SVD sensors 
in a random manner; the standard deviation of these 
shifts and rotations are 15 ^m and 0.15 mrad, respec- 
tively. Effects of tag-side interference [lOl are estimated 
using B — ^- D*£h'. The dominant systematic contribu- 
tions on S are from the uncertainties in the signal and 
background fraction and the resolution function. The 
systematic error on A is dominated by the tag-side inter- 
ference. All these contributions to the systematic errors 
are summed in quadrature to give 0.07, 0.05 and 0.29 for 
S and 0.06, 0.05 and 0.13 for A for the fuU il/^^o^o, MRl 
and non-MRl samples, respectively. 

Ensemble tests are carried out with MC pseudo- 
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experiments using the values of S and A obtained by 
the fit as the input parameters. From 10,000 pseudo- 
experiments, we find that the statistical errors obtained 
in our measurement are within expectations. 

In summary, we have performed a measure- 
ment of the time-dependent CP asymmetry in 
the decay B'^ K^tt'^j with i^^Tr" invariant 

mass up to 1.8GeV/c^, based on a sample of 
535 X 10^ BB pairs. We obtain CP-violation pa- 
rameters SxOjr«-y = — 0.10 ± 0.31(stat) ± 0.07(syst) 
and Ako^o^ = -0.20 ± 0.20(stat) ± 0.06(syst) 
for the full Ks^'^ invariant mass region, and 
5^*0^ = -0.32+°;j^^(stat) ±0.05(syst) and 
Ak'O^ = -0.20 ± 0.24(stat) ± 0.05(syst) for the 
mass region around K*'^{892). This measurement 
supersedes our previous measurement With the 

present statistics, we do not find any significant CP 
asymmetry and therefore no indication of new physics 
from right-handed currents. 
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